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Abstract 



1.1 Context and motivations 



This article deals with some stochastic population protocols, motivated 
by theoretical aspects of distributed computing. We modelize the problem 
by a large urn of black and white balls from which at every time unit a 
fixed number of balls are drawn and their colors is changed according 

Qto the number of black balls among them. The limiting behaviour of the 
. composition of the urn when both the time and the number of balls tend 

to infinity is investigated and the proportion of black balls is shown to 
I ^ I converge to an algebraic number. We prove also that, surprisingly enough, 

not every algebraic number can be "computed" this way. 

> 

m 1 Introduction 

in 

The aim of this article is to tackle some questions of distributed computing 
in theoretical computer science, from a statistical mechanics standpoint. Dis- 
tributed computing deals with large computing systems using many small pro- 
l cessing elements. These small elements are thought as elementary objects in 

^ a complex network whose interactions at a low level may be pretty difficult to 

k>l understand and modelize. There is a clear analogy with statistical mechan- 

I ics, in which physical systems are well described at a macroscopic level, while 

molecular-level phenomena seem chaotic. 

This work is motivated by recent studies in population protocols (see |AR07) 
for a detailed introduction), which are models of decentralized networks con- 
sisting of mobile agents interacting in pairs. The way agents interact is known 
(and assumed to be simple) but not their movements. These movements are 
driven by an "adversary" , which picks, at each time step, two agents according 
to a process only assumed to be fair (roughly speaking, the fairness condition 
ensures that any possible configuration is eventually attained ; see again |AR07) 
for a formal definition). 

Let us be more precise. We are given a finite set S of states, a transition rule 
^ : 5^ — )■ 5, and n > 2 identical agents, which may be at any moment in one of 
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the card(S') possible states. A population protocol associated to ^ is a dynamical 
system {at)tefi on 5*" where at each time step two agents are chosen, and their 
states updated. Updating is made according to (j): if x, y, respectively in states 
e, /, are chosen, then both their state is turned into (j){e, /). Population protocols 
are usually designed to compute predicates over a set of boolean variables. A 
different but related question is addressed in this article: we use population 
protocols to compute real numbers instead of predicates. 



To compute numbers with population protocols, we rely on the classical formal- 
ism of stochastic urn models, where each ball stands for an agent. 

Let us describe more precisely the model we deal with. We fix an integer 
k > 1, a real number xq G [0,1] and a rule f : {0, . . . , fc} — >■ {black, white}. 
At time 0, there are n balls in the urn and each of them is randomly colored 
in black with probability xq or in white with probability 1 — xq, independently 
from the n — 1 others. The triplet (fc, /, xq) is referred to as the rule of the urn. 

At each time unit, k balls are picked randomly, uniformly and independently 
from the past. Let i be the number of black balls among them, then all the fc 
balls are recolored in the color f{i) and put back in the urn. This model is 
a generalization of the famous Ehrenfest urn model (see for instance |Dur91) ^ 
which corresponds to fc = 1 and i— ^ black 1 i— > white (or equivalently at each 
time step a single ball is picked and its color is changed). 

Definition 1. The real number a is said to be computable if there exists a rule 
(fc, /, Xq) such that, for any e > 0, there exists c > such that 



where is the proportion of black balls at time £, and [fcj stands, as usual, 

for the largest integer smaller than fc. 

Roughly speaking, it means that if we start with a large number n of balls 
and wait for a linear time in n (larger than cn), the proportion of black balls is 
with high probability close to a. 

Before stating our main results, let us first discuss some important features 
of our model. 

• The results we aim for are about the asymptotic properties of population 
protocols, when the size of the population grows to infinity. From this 
perspective, the choice of picking agents uniformly at random appears to 
be the natural generalization of fairness to large (or even infinite) popu- 
lations. 

• It is important to point out the strong assumption that the fc balls are all 
turned into the same color. From our original setting, this is motivated by 
the fact that in the complex network there is no hierarchy and not much 
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communication between the agents: when they meet, they instantly ah 
take the same decision. 

• Because the number we compute corresponds to the proportion of black 
balls, it is sufficient to consider only two possible states for the balls. This 
is why we focus on this case. 

1.3 Description of the results 

The results of this paper are twofold. We first describe in Section [2] the asymp- 
totic behavior of the proportion of black balls for any fixed rule. In Proposi- 
tionjsj the process {Xy^^^ )t>o is shown to be close to the solution of an ordinary 
differential equation. We then characterize in Theorem [4] the numbers com- 
putable by a given rule as the roots of a polynomial related to this differential 
equation. 

Section [3] is devoted to studying the set of computable numbers. This relies 
on a combinatorial analysis of the latter polynomial. In Theorem [5j we show 
that the set of computable numbers is dense in [0,1]. On the other hand, 
we prove in Proposition [7] that, surprisingly enough, rational numbers are not 
computable (except for a few explicitly exhibited ones). 

We conclude this introduction by comparing our results with related works. 
In 'BCC+09I, the black/white case with k = 2 has already been handled but 
with differences in the approach and the statements of the results. The main 
difference is that the authors of |BCC"'"09] had to assume that the initial pro- 
portion xq of black balls is close to the number a one wants to compute]^ (as a 
counterpart, the main result of BC C"'"09] gives a interesting and precise descrip- 
tion of the fluctuations of ^fe"'' around its mean, for k large). The techniques 
developed in the present paper allow us to free ourselves from this restrictive 
assumption. 

Our main results (Theorem [5] and Proposition [?]) may seem surprising as we 
might expect that any algebraic number would be computable in our setting. 
This should be put in perspective with a simultaneous and interesting result by 
Bournez, Fraigniaud and Koegler |BFK12j . They study stochastic population 
protocols in which agents are only picked in pairs but the number of possible 
colors is arbitrarily large and the new colors of the two balls may be different. 
In this different setting, any algebraic number is computable. Note that, as in 
|BCC"'"09l , the results in [BFK12j hold only with the assumption that the initial 
proportion of black balls is close to an equilibrium. 

Let us also note that the link between the evolution of some stochastic pop- 
ulation protocols and that of an associated ordinary differential equation has 
been used for the first time by Chatzigiannakis and Spirakis [CS08. in a some- 
what different context; they study some qualitative properties of the differential 
equation in order to discuss the stability of the underlying protocol. 

1 Indeed, the hypothesis - h{x)\ 0" in f |BCC+09) . Th.2) imphes that "1 - 

2(Xq"' )^ + Xg"^ ^^-^ — > 0" and thus the initial proportion has to go to a. 
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Before going through details, we present now a simple introductory example. 



1.4 Heuristic : the example of (3 — v^)/2 

Take fc = 2 and consider the function / : i— > black ; 1 n- white; 2 !—> white, as 
illustrated below 




Recall that X^^^ denotes the proportion of black balls at time £. The sequence 
{X^^^ ) defines a Markov chain on the set ■ ■ ■ , ^} , which admits a unique 

invariant measure Tr'"^ Transition probabilities of the chain X are clearly ra- 
tional numbers so the components of 7r("\ as a solution of a linear system of 
rational equations, are rational numbers. Its mean is thus rational. Hence 
ergodic theorem for Markov chains states that almost surely: 



in) 



X 



in) 



in) 



Mean(7r(")) e 



To get a hint for the asymptotic behavior of p^^"^ , let us compute the conditional 
expectation of the increments of xl"-* : 



E 



n 
2, 



' (both balls are white) 
' (both balls are black) , 



+ 



2 ("7") 



1 nx{n 
n f! 



' (one ball is white, one is black) 



/nx\ 
\ 2 ) 



{2{l ~ x)^ - 2x{l - x) - 2x^) 



(1) 



Take now £ large, our system converges to its stationary regime, and thus we 
expect the righ-hand term in ([!]) to vanish. Hence, for large n, p*^"' should be 
close to the irrational number (3 — \/5)/2 0.382 . . . , which is the only root 
of the polynomial 2(1 - X)"^ - 2X{1 - X) - 2X^ in [0, 1]. We let the balls 
"compute" (3 - %/5)/2. 



2 Limiting behavior of urns 

We study in this section the sequence X^^"^ :— (x|"-')£>o of the proportions of 
black balls in the urn. As mentioned above, this sequence is a Markov chain 
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with state space ^, ■ ■ ■ ,^}- We denote by Ef the set 

Ef ^{0<i<k; f{i) = black} , 

recall that the triplet (fc, /, xq) (or equivalently, the triplet (k, Ef^xo)) is the 
rule of the urn. 

Following the heuristic of Section 



1.4 



we associate to the rule (fc, /) the 



polynomial b — bf defined by 



b{y) 



f 

E 



leEf 



ky\l-yf-'-ky 



\k — i 



= kV{Bk,yeEf)-ky, 

where Bk,y is a binomial random variable with parameters {k,y). In the above 
example (fc = 2 and Ej = {0}), this gives 

2P(B2,, =0)-2t/-2(l-2/)2-22/, 

which is indeed equal to the polynomial in ([T]). 

The meaning of b{y) can be understood as follows: when n is large, picking 
fc balls uniformly among n balls, a proportion y of them being black, almost 
amounts to perform fc times an experiment with a probability y of success. The 
quantity b{y) then represents the expectation of the evolution of the number of 
black balls, after putting back the fc recolored ones, as stated in the following 
lemma: 

Lemma 2. For y ^ c/n with c E {0,1,2,..., n}, set 

b(-\y)^E[x[-^-X^-^\xt^^y 

The map y nfo*-"-* (y) converges uniformly to b{y) on the set {0, 1/n, 2/n, . . . ,n/n}. 
More precisely, for n big enough, 



max 

0<c<n 



n6(")(c/n) - b{c/n) 



< 5¥ 



Proof. The probability that i balls are black when fc balls are picked in a urn 
that contains n balls among which ny are black is equal to: i^^^f) (^) , 
indeed to pick exactly i black balls, we need to pick i black balls among the ny 
black balls and n—i white balls among the n—ny white balls. For y € {0/n, n/n} 
there is nothing to prove since n&'^"-*(0) — 6(0) ~ 7i6'^"-'(l) — 6(1) = 0. For any 
y S {^/n, 2/n, . . . , (n — l-)/n} we write 



n6(")(y)-6(y) 



fc 

E (fci'SB/ 

k 

E 

1=0 



(ny\ fn-ny 
\ i J\ k-i 

(fc) 



- y) 




y'ii-y) 



k—i 



(T)(t: 

)Qy'i^-y)' 



k — i 



- 1 
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Let us handle the last term 



^ (7)(T-?) _^ (T)»! (T-7)(fc-0! 

To show that ^ goes to zero (and hence, so does nb'^^\y) — h{y)), recall that, 
when j is fixed and m goes to infinity, (™) ~ e'^mP / j\ : the last three terms 
in ([2]) converge to one. 

In order to prove that the convergence is uniform, a little more work is 
needed. First, recall that a consequence of the Stirling formula is that, for any 
integers j < m, 

2p ("'^ 



exp [-^)< 



< 1. 



(3) 



Plugging this in ^ gives that 

(nv\ (n-ny 
\ i J \ k-i 



1 - exp(-2fcVri) < 1 
and thus 

1 - 



— < 1 — expi ) exp( 

1){l)y'{l-y)^-' ~ ny' n-ny ' 



rny\ /n-ny 

\ 1 / \ ~ 1 



< 1 — exp I — 



If ^1 \/n < y < I — l/\/n, then this last quantity is less than 4k'^/y/n. Then, for 
any y G [1/ ^/n; 1 — 1/ y/n\ , one has 



- b{y)\ < J2 - z| . W'(l - y)^ 



nv\ (n-ny\ 
i J \ k-i I 



We now deal with the case where y < and prove that both nb''"'\y) and 

h{y) are close to fcloe£;^(l — y)^ ■ We assume here that n > 4fc^ and therefore 
y < l/2k . 



nb^''Hy)-kloeEf{l-y)'' = klo 



deEf 



ny\ /n-ny\ 
y V fc-O / 



4=1 



■iy\ /n-ny\ 
i ) \ k-i ) 

it) 



Then, using again ([s]) and inequalities 1 — fcy < (1 — y)^ < 1 — ^ and 1 — a; < 
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e-^ < 1- f for 0<a;< 1/2, 
|n6(")(y) - kloeEfil - < kloeEf 

< kloeEf 



/ n — ny\ 

\ k J (A ^.\k 



/n-ny\ 

^ fc ^ - fl - v)'' 



fe (ny\ (n-ny\ 
l,Sr^ V i A k-i J 

i=l \k) 



< fc(l - exp(-2fcV(" - - 1| + ^^(1 - (1 - j/)'=e-2^'/("-"^') 

41,2 

<fcxlx — + k(l-(l- ky){l - Ak'^/n)) 

n 



< Sk'^/n + k-'y 



2^ Ayk'^ ^ 5fc2 
n ~ \/n 



On the other hand, 



\h{y) - fcloeE.Cl - = I E (^I'ss/ ' • - 

< A;(l - (1 - yf) < k^y < k'^/s/n. 
This proves that for any y < 1/ y/n, 

|n6("' (j/) - b{y) \ < 5k'^/s/n + k'^/Vn < 6A;V^/n, 
and concludes the case y < 1/ y/n. The case y > 1 — 1/ \fn is symmetric. □ 

The rest of the section is devoted to the study of the convergence of X("\ 
For that purpose, wc define t n- x{t) as the unique maximal solution of the 
ordinary differential equation (ODE) such that 

x' = h{x) and x(0) = xq (4) 

(recall that xq is the initial proportion of black balls in the urn). First notice 

that since 

6(0) = fcP (Bfc Q^Ef)= klo^Ef > and 6(1) = A;P {Bk kk € Ef)-k = kikes f-k < 0, 

(5) 

this maximal solution Xt actually remains in the interval [0, 1]. To describe; the 
asymptotic behavior of the sequence X^"', we speed it up by a factor n, by 
setting Xn{t) = and prove that Xn{t) is well approximated by x{t) when 

n is big: 

Proposition 3. For each rule and for any real numbers to,E > 0, there exist 
A,B>0 such that, for each n> 1, 

P (sup \xn{t) - x{t)\ >e] < Ae-^". 

\t<to J 
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Proof. As this proposition can be seen as an instance of the general theory of 
large deviations for Markov processes, sometimes known as Kurtz's Theorem 
(see |SW95j ). we only outline the main ideas of the proof. 

(n) 

For sake of conciseness we set Xk := Xj^ and introduce the classical mar- 
tingale 

k-l 



Mk=Xk-X^-Y,U^'>{X,). 



1=0 



This equation enables to rewrite X\j^t\ as: 



" [_nt\ I n 

XL„tJ =Xo + Myr,t\ + / nh^'^'^ {Xyns\ )ds. 



Then 



Xy^t\ -Xt = M L„tj + / nM") (Xl„,j ) - &(Xl„,j )ds 



\_nt\ /n 

b{X[nsi) - b{xs)ds + {Xq - Xo) + (a^Lntj/n - Xt), 

since a; L„tj /„ — xq+J^^*^^^ b{xs)ds. Now, in order to bound, /(t) := sups<( |a;„(s) — 
x(s)|, we write 

p \_nt] / n j* [nt\ j n 

fit) < sup |Ml„,j 1+ / \nb^"\x^{s))-b{x,,{s))\ds+ / \b{xn{s))-b{x{s))\ds 

s<t JO JO 

+ (Xo - Xq) + sup \xins\/n ~ Xs\. 
s<t 

The probability for the first term to be large can be bounded with a concen- 
tration inequality for martingales while the second term is bounded thanks to 
Lemma [2] The probability that the two last terms are greater than e is as small 
as desired. Since the third term is smaller than sup J /(s)c?s, an applica- 
tion of the Gronwall Lemma gives a bound for V{f{t) > e). Again, we refer to 
f |SW95j . p.76-84) or f |DN08j . p.45-46) for details. □ 

The Cauchy-Lipschitz Theorem implies that x'(t) is never equal to zero 
(unless X is constant and equal to a root of b). Hence any solution of Q is 
monotonous, and converges to a root of b. If b{xo) > (resp. < 0) then the 
solution starting from xq converges to the smallest (resp. largest) root of b 
greater (resp. smaller) than xq, denoted by a. We gather this observation with 
Proposition [3] to obtain our main result: 

Theorem 4. Assume that b{xQ) > (resp. < 0) and let a be the smallest (resp. 
largest) root of b no smaller (resp. no greater) than xq. For any e > 0. there 
exist some constants c > and A, B > such that for each n 



x\''\ 

[cnj 



> e < Ae 
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In particular it implies that the rule {k, Ej , xq) computes the number a, as 
defined in Definition^ 

As regards a more quantitative aspect on time and space complexity, we 
point out the recent article ABll in which this question is discussed for the 
case of fc = 2 and a = 1/ \/2 (but the method extends to other situations) . 

Proof of Theorem First note that such an a always exists by ([s]) . Take now 
c large enough, so that \x{c) — a\ < e/2. It suffices then to write 



lcn\ 



>e 



)<v{\xl:l^-x{c)\>e/2) 



Proposition [3] gives the desired bound for the right-hand side. 



□ 
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Figure 1: Two simulations of with n = 2000 balls up to time 2000, with 

the flow of the corresponding ODE. Here, fc = 8, £' = {0, 4, 6, 8} and xq = 0.5. 
The polynomial h is equal to 8 ((1 - xf + 70x^(1 - xf + 28x^(1 - x)'^ + x^) - 
8a; and the corresponding a is approximately equal to 0.2079. 



3 The set of computable numbers 

We give in this section some properties about the set C of numbers that can be 
computed by our urns. A first basic observation is that each element of C is the 
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root of a polynomial 6/ and hence is algebraic. Moreover we have the following 
properties: 

Theorem 5. The set L 

(i) is symmetric with respect to 1/2 ; 

(ii) is dense in [0, 1] ; 

(iii) contains numbers of any algebraic degree ; 

(iv) does not contain every algebraic number. 

Proof, (i) Let a in £ and {k,E,xo) a be rule computing to a. We denote by 
E* the set defined by: 

i e E* ^ k - i ^ E, . 
Let b* be the polynomial associated to the new rule {k,E*, 1 — xq), we have 

b*{l-a)^kF{Bk,i-aeE*)-k{l-a) (6) 
^ -k¥ (i3fe,„ eE) + ka. (7) 

Hence, 1 — a is a root of b* . One checks easily that if the solution of the ODE 
y' = b{y),y{0) = yo converges to a, then the solution oiy' = b*{y), y{0) = 1 — yo 
converges to I — a. 

(ii) Let a/b be a rational number in [0, 1], and e,d two positive reals such that 

{a/b - e, a/b + e) C {6,1- S). 

We are looking for a number a £ {a/b ± e) and a rule {k, E, xq) such that the 
associated ODE converges to a. In particular it is necessary that: 

F{Bk,c. eE)^a. (8) 

Fix for now the integer fc, and consider the set 

Ea.b = {i < fc; i = 0, 1, 2, . . . , o - 1 (mod b) } . 

The proof relies on the following lemma: 

Lemma 6. For any < 6 < 1/2, there exists A > such that for any integer k 
and X e {6,1 — 6), 

|P {Bk^^ = 0, 1, . . . , a - 1 (mod b) ) - a/b\ < e'^K (9) 

Proof of Lemma A proof based on linear algebra would give the best constant 
A. As we do not need here this exact value, we give a probabilistic and shorter 
proof. The value modulo {b — 1) of a random variable Bk.x is the position at 
time k of the walk X = {Xi)fyo on {0, 1, . . . ,b — 1} starting from Xq — and 
with probability transitions 

V{Xi,+i ^Xt + 1 (mod b)) = l- f{Xi+i = Xi (mod b))=x. 
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Figure 2: A plot of the maps x ^ x and x i-^ f'{Bk,x G Ea.b), for k ~ 30, 
a/6 = 3/5. 



starting from Xq — 0. It is clear that this Markov chain admits as unique 
stationary measure the uniform measure tt over {0, 1, . . . , 6 — 1}. By the general 
coupling inequality (see jLin92| Chap. 1. 2.), the desired quantity is smaller than 

where X, X are two i.i.d. copies of X, starting from and tt. These two walks 
meet necessarily if during b successive steps X goes b steps forward while X 
remains motionless. This occurs with probability x^{l — x)'', hence 

|P (6fc,, = 0, 1, . . . , a - 1 (mod b) ) - a/b\ < (l - x'il - x)') ^'^'^ , 

which decays exponentially in fc, provided x is bounded away from and 1. □ 

Assume fc is a multiple of b, this ensures that G E and 1 ^ E and thus that 
neither or 1 is a root of b. So we might as well take a smaller 6 such that all 
the roots of b in the interval [0, 1] belong in fact to {S, 1 — 5). Let k be such that 
g-Afc ^ ^ _ Q g_ rp-j^g solution of y' — b{y) starting from 0.5 converges to 

a root of b. By Lemma [6] such a solution belongs to (a/6 ± e) (see Figure [2]). 

(iii) Fix fc > 1 and consider the set E = {1}. The associated polynomial is 

fca(l - a)''"^ - a. 

Its unique root in (0, 1) is 

xa = 1- ''^^/l/k, 
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which has algebraic degree fc — 1 . 



(iv) We give in fact in the next proposition a much stronger result stating that 
almost no rational numbers belong to the set C 

□ 

Proposition 7. Let x = p/q be a rational number such that gcd(p, g) = 1 and 
g > 4 then x (Ji C. 

Before proving this proposition, observe that the only rational numbers be- 
tween and 1 that do not satisfy the above conditions are 0, 1, 1/2, 1/3 and 
2/3. These numbers all belong to C and are respectively computed by the rules 
(1,0,0.5), (1, {0,1}, 0.5), (2,{1},0.5), (3, {0,3}, 0.5) and (3, {1, 2} , 0.5). 

We proceed by contradiction. Let x = p/q such that gcd(p, q) = 1 and g > 4 
and assume that p/q G C. Since it implies that 1 — p/q G £, we can assume 
without loss of generality that p > 3. Let {k,E,xo) be one of the rules that 
admits p/q as a solution, we can hence write 

J2(^)p'(i~p)'''-pi'''- (10) 



We now use some well-chosen reductions modulo p to deduce from this 
relation that fc = 1 (mod p") for every n, which leads to a contradiction 



(take for example n equal to fc). Reduction of (10) modulo p implies that 
lo^Eq'' = (mod p) , which yields ^ E, since gcd(p, q) = 1. We go one step 
further, reducing (|10[) modulo p^ and dividing by p leads to the relation : 



lieEHq - pf ^ =q'' ^ (mod p) , 

in which the left-hand side can be simplified into ligsfcg'^^^ (mod p) . Since 
gcd(p, q) = l,we obtain lig^jfc = 1 (mod p) , from which we readily deduce that 
1 € E and fc = 1 (mod p) . 

We now proceed by induction to show that fc = 1 (mod p") , for every n < k. 
The following lemma will be useful: 

Lemma 8. Assume fc = 1 (mod p"^^) , with 2 < n < k. Then for any 2 < i < 
n, (*^) = (mod p"^*+i) . Moreover if p ^ 2 (mod 4) , the result is also true 
for i = 2. 

Proof. It is enough to prove the lemma for p being a power of a prime number, 
otherwise writing the decomposition of p = 11 -Pf' i'^^o a product of prime 
numbers and applying the result for each of the term gives the result. We 
start with the classical relation: 

- 2 



i(i- 1)[^, ) = fc(fc- 1) 
Since fc = 1 (mod p"~^) , we get 



gcd(i,p""i)gcd(i- l,p"-i)( ^ ) =0(modp"-i). (11) 
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Now, for i > 4, we have that i < 2'~^ < and hence gcd{i,p"~^) 

gcd(i,p*^^), smce we assume that p is a power of a prime. Similarly gcd(i 
— gcd(i — l,p'~^). Now, since gcd(i,i — 1) = 1, 



Then ( 11 ) yields the desired resuh. The cases i — 2 and i = 3 are dealt with a 
direct computation. □ 

To continue the proof of Proposition [7j we need to proceed differently de- 
pending on p being or not equal to 2 modulo 4. Assume first that p ^ 2 (mod 4) 
and that we proved k = 1 (mod p"~^) for some n > 2. To carry on the recursion, 
we write the reduction of (10) modulo p""*"^, 

kpiq - p)^-^ + UeE (fjp\q - P)''-' ^ pq'-' (mod p^+') . 

i=2 

Lemma |8] implies that each term in the sum of the l.h.s. vanishes. Expanding 
the remaining term kp{q — p)^~^ and dividing both sides by p gives 

fe-i /, 
k 



i=l 



(1 - k)q^-^ + Y^{k ]p'q'"~^~' = (mod p") 

(1 - k)q''-^ EE (mod p") 



This proves by induction that k = 1 (mod p") for any n, which leads to a 
contradiction and concludes the proof in this case. 

Assume now that p = 2 (mod 4) and fc = 1 (mod p"^^) for some n > 2. We 
first observe that 



q' 



Thus ( 10 ) can be written as 

E(-)p^(9-^')'"' = 9'-m'"'- (12) 

Taking the reduction of the latter equation modulo p"+^ and using that (^)p* = 
(mod p^^^) when i > 2 gives 

{q - pf + l2iE Q)p'(<Z - P?-^ ^q'- pq'-' (mod p"+l) 

q' - kq'^-'p + ( 'f] q'-V + l2^E ( %'iq pf'^ ^ q" - pq''' (mod p"+l) 



2/" " ^^"'V2 
k 
2 



(1 - fc)pg'=-i + Cf) q'^-V + Me - Vf-'' ^ (mod p"+i) . (13) 



13 



We focus on the last term of the left-hand side and write: 



^l2^E(^yV-Hmodp-+'), 



since p/2 e N and (fc - 1)^^ = (mod 
Then ( 13 1 divided by p can be written 



(1 + (l + l2^^,)Q) g'-'p- 



= (mod p") 



If 2 ^ E, the second term disappears and we are left with (1 — k)q'^ ^ = 
(mod p") and hence fc = 1 (mod p") , otherwise we get 

Y)^0(mod p"). 

We conclude the proof by noticing that (<? + ^) is both prime with p/2 and 
odd, hence prime with p. 



4 Conclusion 

These first results raise some interesting theoretical questions in the current 
research on the computational power of population protocols. Although our 
model is very general and allows to compute a large set of numbers, some 
algebraic numbers as "simple" (on a computational point of view) as 1/5 are not 
computable. A natural question is then to ask if the set C has a nice structure 
: has it interesting symmetries? can it be endowed with a certain algebraic 
structure which is consistent with computability? In other words: does there 
exist an operation (g) such that if x and y belongs to £ then a certain combination 
of their associated rules computes x® yl 

As already mentioned, it is proved in |BFK12| that any algebraic number is 
computable for k — 2 and q > 2 colors, but with the significant difference that 
the 2 balls may be turned into two different colors. A question remains: what 
happens in our model with k > 2 and two states if we consider more general 
rules for which the k balls may be recolored differently from each other? With 
this new model it is possible to compute any rational number but we still do 
not know if any algebraic number is computable ; [BFK12j suggests that this 
should be the case. It also would be interesting to study whether adding more 
colors but still requiring that the k balls all turn into the same color has a bigger 
computational power. 
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